o 

U 

Or 
(N 

w 

X 

43 
Oh: 

6 

CO 



> 
in 

(N 

o 



3 



Pulse phase dependent variations of the cyclotron absorption 
features of the accreting pulsars A 0535+26, XTE J1946+274 and 

4U 1907+09 with Suzaku 

Chandreyee Maitra 1 and Biswajit Paul 
Raman Research Institute, Sadashivnagar, Bangalore-560080, India 

cmaitraOrri . res . in ; bpaul@rri . res . in 
Received ; accepted 



1 Joint Astronomy Programme, Indian Institute of Science, Bangalore-560012, India 



-2- 



ABSTRACT 

We have performed a detailed pulse phase resolved spectral analysis of 
the cyclotron resonant scattering features (CRSF) of the two Be/X-ray pul- 
sars A0535+26 and XTE J1946+274 and the wind accreting HMXB pulsar 4U 
1907+09 using Suzaku observations. The CRSF parameters vary strongly over 
the pulse phase and can be used to map the magnetic field and a possible devi- 
ation form the dipole geometry in these sources. It also reflects the conditions 
at the accretion column and the local environment over the changing viewing 
angles. The pattern of variation with pulse phase are obtained with more than 
one continuum spectral models for each source, all of which give consistent re- 
sults. Care is also taken to perform the analysis over a stretch of data having 
constant spectral characteristics and luminosity to ensure that the results reflect 
the variations due to the changing viewing angle alone. For A0535+26 and XTE 
J1946+274 which show energy dependent dips in their pulse profiles, a partial 
covering absorber is added in the continuum spectral models to take into ac- 
count an additional absorption at those phases by the accretion stream/column 
blocking our line of sight. 

Subject headings: X-rays: binaries- X-rays: individual: A0535+26- individual: XTE 
J1946+274- individual: 4U 1907+09- stars: pulsars: general 
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Introduction 



The X-ray binary sources in which the compact object is a highly magnetized neutron 
star, often with a massive companion are called accretion powered pulsars. Due to the strong 
magnetic field of the neutron star, the matter here flows along the magnetic field lines to the 



poles of the system, forming an X-ray emitting acc retion column above it ( jPringle fc Rees 



1972 



Davidson fc Qstriker 



1973 ; 



Lamb et al. 



19731 ). Another important consequence of the 



strong magnetic fields (~ 10 12 G) are the cyclotron resonant scattering features (CRSFs) 
formed by the resonant scattering of photons by the electrons which are quantized into 



Landau levels forming absorption like features at multiples of E c = 11.6keV x yp 



x 



B 



10 12 G' 



E c being the centroid energy, z the gravitational redshift and B the magnetic field strength 
of the neutron star. The CRSFs thus provide a direct tool to measure the magnetic 



field strength o 



(jTrumper et al 



the neutron star. It was fi rst discovered in the spectrum of Her X-l 



1977 



discovered so far (jPottschmidt et al. 



Truemper et a 



1978) and about 20 sources with CRSFs have been 



20121 ). The CRSFs which are found mostly in high 



mass X-ray binaries, about a half of which are transient sources, lie between the energy 
range of 10-60 keV. In addition to the magnetic field strength, the CRSFs also provide 
crucial information on the emission geometry and its physical parameters like the electron 
temperature, optical depth etc. 

Pulse phase resolved spectroscopy of the cyclotron parameters is an especially useful tool to 
probe the emission geometry at different viewing angle as the neutron star rotates. It can 
further be used to map the magnetic field geometry of the neutron star. Since the CRSFs 
also show variations with luminosity and spectral changes, to perform pulse phase resolved 
analysis, care should be taken to obtain the results solely due to the changing viewing 
angle by averaging over the data stretch with similar counts and spectral ratios. Proper 
continuum modeling of the energy spectrum also plays an important role in phase resolved 
analysis. Suzaku, with its broadband energy coverage is most ideally suited in this regard. 
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A0535+26 is a Be /X-ray binary pulsar which was discovered during a giant outburst 



in 1975 by Ariel V ( [Rosenberg et al 



19751 ). It consists of a 103 s pu lsating neutron star 



is ~ 2 kpc (jGiangrande et al 



1980: 



with a 09.7IIIe optical companion HDE24 5770 ( 

of e=0.47, with orbital period of 11 1 days (jFinger et al. 



Bartolini et al 



Steele et al. 



19781 ) in an eccentric orbit 



20061 ) . The distance to the source 



19981 ). Upto six giant outbursts have 



been d etected in this source so far, the latest ones during 2009/2010 (iCaballero et al. 



2011a|JbJ). The last giant outburst was followed by two smaller outbursts with a periodicity 
of 115 days which is longer than its orbit al period. Precurso rs to the giant outburst was 



also observed with the same periodicity (IMihara et al. 



20101 ). CRSFs at ~ 45 keV and 



100 keV were discov ered in this source during the 1989 giant outburst with HEXE 



( iKendziorra et al. 



19941 ) . 



during the 1994 outburst (IGrove et al. 



he second harmonic at ~ 110 keV was confirmed with OSSE 



19951 ). although the presence of the fundamental 



at ~ 45 keV was dubious. It was later confirmed durin g the 2005 outburst with Integral, 



RXTE ( ICaballero et al.l 120071 ) and Suzaku observations ( ITerada et al 



20061 ) 



XT 



J1946+274 is a transient Be/X-ray binary pul sar discovered by A SM onboard 



RXTE (ISmith et al. 



19981 ) . and CGRO onboard BATSE (IWilson et al. 



19981 ) during a giant 



outburst in 1998, revealing 15.8 s pulsations. The optical counter part was identified as a n 



optically faint B ~ 18.6 mag, bright infrared (H ~ 12.1) Be star (jVerrecchia et al.l 120021 ). 



'he source 



( IPaul et al- 



ias a moderately eccentr ic orbit of 0.33 with an orbital period of 169.2 days 



2001 



Wilson et al 



20031 ). After the initial giant outburst and several short 



outbursts at periodic in tervals, the source went into quiescence for a long time until the 



recent outburst in 2010 (ICaballero et al 



2010b|). A CRSF was discove red at ~ 35 keV from 



the RXTE data of the 1998 outburst observations ( iHeindl et al. 



200lh . 



4U 1907+09 is a persistent wind 



in the Uhuru surveys ( IGiacconi et al 



accre t ing high mass X-ray binary discovered 



1971 : 



Schwartz et al. 



19721 ). having a highly 



redenned companion sta r (08-09 



M = 7 * 1Q- 6 M yr~ l ( 



Cox et a 



8.3753 days (lin 't Zand et al 
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a) of magnitude 16.37 mag and a mass loss rate of 



20051 ). It has a moderately eccentric (e=0.28) orbit of 



1998). It is a slowly rotating neutron star with period of ~ 440 



s, and has showed several episodes of torque rever sals with a steady spin-down from 1983 to 



1998 ( 



in 't Zand et 



2003 ( Bavkal et al. 



al 



1998 



Mukerjee et al. 



20011 1 . a much slower s pin- down from 19 98 to 



20061 ). a torque reversal bet ween 2004 and 20 05 (jFritz et al.ll2006l ) and a 



second torque reversal between 2007 and 2008 (jlnam et al. 



20091 ). which restored the source 



to the same spin-down rate be 



from the Ginga /observations (JMakishima &: Mihara 



'ore 1998. A CRSF at ~ 1 9 keV was reported usin g data 



1992 



Makishima et al 



19991 ) and was 



later con firmed from the Beyyq S aX observati o ns wi th the discovery of a harmonic at 



36 keV. (ICusumano et al. 



1998 ). 



Rivers et al. 



(120 101 ) performed a time and phase resolved 



analysis of the Suzaku observations of the source made during 2006 and 2007. 
Here we present the results obtained from a pulse phase resolved spectroscopic analysis of 
these three sources with a motivation to investigate the variation pattern of the cyclotron 
parameters with pulse phase. The pulse phase dependence of the CRSF parameters are 
presented for the first time for A0535+26 and XTE J1946+274, whereas a more detailed 
result is pre s ented for 4U 1907+09 which is in agreement with the earlier results of 



Rivers et al. 



LO). 



(|2010| ). The analysis is done taking into account various factors which might 



smear the pulse phase dependence results. The results presented in this paper are one of the 
most detailed results on pulse phase resolved measurements of CRSF available so far. 

2. Observations & Data Reduction 



There are two sets of scientific ins truments onboard Suzaku. The X-ray Imaging 



Spectrometer XIS (jKoyama et al 



20071 ) consisting of three front illuminated CCD detectors 



(FI :XIS0, XIS2, XIS3) and one back illuminated CCD detector (BI: XIS1) work in the 
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0.2-1 2 keV range and the Hard X-ray Detector (HXD) made of PIN diodes (ITakahashi et al. 



20071 ) and GSO crystal scintillator detectors cover the energy bands of 10-70 keV and 



70-600 keV respectively 



Suzaku ( JMitsuda et al. 



120071 ) observed A0535+26 twice, once on September 14-15 2005 
during the decline of the second normal outburst of 2005, and again on August 24 2009 
during the decline of the 2009 normal outburst. We have chosen the 2009 observation (Obs. 
Id-404054010) for our analysis because of the longer duration (exposure ~ 52 ks) and its 
'HXD nominal' pointing position which is more suitable for CRSF studies, although the 
count rates were comparable for both the observations. The XIS's were operated in the \ 
window' 'burst' clock data mode which has a total time resolution of 2 s. 
XTE J1946+274 was observed on October 11 2010 (Obs. Id-405041010) just after the 
peak of the September/October 2010 normal outburst. The source was observed for ~ 51 
ks in the 'HXD nominal' pointing position, and the XIS's were operated in the '| window' 
'normal' clock data mode which has a time resolution of 2 s. 

4U 1907+09 was also observed twice with Suzaku, once on May 2006, and again on April 
2007. We have chosen the 2007 observation (Obs. Id-402067010) for our analysis because 
of similar reasons as in the case of A0535+26, i.e. longer exposure of ~ 158 ks and the 
'HXD nominal' pointing position. The XIS's were operated in 'normal' clock data mode 
with no window option which has a time resolution of 8 s. The XIS data was reduced 
and extracted from the unfiltered XIS events, which were reprocessed with the CALDB 
version 20120428. We checked for any significant photon pile-up effect in the reprocessed 
XIS event files. To perform pile-up estimation, we examined the Point Spread Function 
(PSF) of the XISs and obtained the count rate at the image peak per CCD exposure as 
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given by Yamada & Takahashi ( lYamada et al.ll201 2rl. Crab data is assumed to be free from 



pile-up and has a value of 36 ct/sq arcmin/ s/CCD exposure at the image peak. Following 
their procedure, the XIS data of A0535 +26 and XTE J1946+274 had values of 2-3 
ct/sq arcmin/ s/CCD exposure at the image peaks, and showed no evidence of significant 
pile-up. 4U 1907+09 on the other hand showed a case of moderate photon pile-up. The 
value obtained at the image center was higher than the Crab Nebula count rate of 36 
ct/sq arcmin/ s/CCD exposure. The radius at which this value equals 36 in the PSF is 
about 15-16 arcsec, and hence 15 pixels were removed from the image center to account 
for this effect. For the extraction of XIS light curves and spectra from the reprocessed XIS 
data, a 4 diameter circular region was selected around the source centroid for A0535+26 
and XTE J1946 +274, and an additional central 15 (~ 16 ) pixels were removed in the 
case of 4U 1907+09 to discard the maximum pile up affected regions. Background light 
curves and spectra were also extracted by selecting regions of the same size away from the 
source. The XIS count rate was 3.6 c/s, 3.1 c/s and 8.1 c/s for A0535+26, XTE J1946 
+274 and 4U 1907+09 respectively. 4U 1907+09 had a ~ 12 % loss in count rate after the 
removal of the photons from the central region due to the pile-up correction. Response files 
and effective area files were generated by using the FTOOLS task 'xisresp'. The HXD/PIN 
light curves and spectra were extracted after reprocessing the unfiltered event files q The 
HXD/PIN background was created by adding the simulated 'tuned' non X-ray background 
event files (NXB) corr esponding to the month and year of the respective observations 



Fukazawa et al. 



( 12009 



j to the the cosmic X-ray background, which was simulated as 



1 http://www-utheal.phys.s.u- tokyo.ac.jp/ yuasa/wiki/index.php 
/How_to_check_pile_up_of _Suzaku_XIS_data 

2 http://heasarc. nasa.gov/docs/suzaku/analysis/hxd _repro.html 
3 http://heasarc. nasa.gov/docs/suzaku/analysis/ pinbgd.html 



suggested by the instrument team_| after applying appropriate normalizations for both 

cases. 

facility 



cases. The corresponding response files were obtained from the Suzaku guest observatory 



3. A 0535+26 and XTE J1946+274 

3.1. Timing analysis: Light curves Sc Hardness ratio & pulse period 

determination 

We performed timing analysis after applying barycentric corrections to the event data 
files using the FTOOLS task 'aebarycen'. Light curves were extracted with a time resolution 
of 2 s for the XISs (0.2-12 keV), and 1 s for the HXD/PIN (10-70 keV) respectively. For 
XTE J1946+274 which has a short pulse period, light curve with the resolution of 10 ms 
was extracted from the HXD/PIN data to search for the pulse period. We applied pulse 
folding and \ 2 maximization technique to search for pulsations in the XIS data for A 
0535+26 and PIN/HXD data for XTE J1946+274. The best estimate of the period was 
found to be 103.47 ± 0.09 s for A0535+26. This value is consistent with the pulse period 
determined from the IN TEGRAL IBIS data during the same outburst at MJD 55054.995 



(jCaballero et al 



2010a]) assuming the spin down value determined from the same. For XTE 
J1946+274, the best-fit period was estimated to be 15.75 ±0.11 s. Orbital correction of the 
pulse arrival times was not required for both the sources having a very long orbital period. 
The XIS and PIN light curves of the sources binned with its pulse period for A 0535+26 
and 10 pulsar periods for XTE J1946+274, are shown in Figure [U The light curves show 



4 http://heasarc. nasa.gov/docs/suzaku/analysis/pin_cxb. html 

5 http:// heasarc.nasa.gov/docs/heasarc/caldb/suzaku/ and used for the HXD/PIN spec- 
trum 
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more or less constant count rate, and do not have any particular trend of variation. For 
each figure, the third panel shows the hardness ratios (ratio of PIN counts to XIS counts) 
which is also more or less constant throughout the observation duration and does not have 
any signatures of spectral variability which might affect the results of pulse phase resolved 
spectroscopy. 

3.2. Energy Dependence of the pulse profiles 

We created the energy resolved pulse profiles for the entire stretch of observations by 
folding the light curves in different energy bands with the obtained pulse period. The pulse 
profiles in the energy range of 0.3-12 keV were created using all the three XISs (0, 1 & 3), 
and in the 10-70 keV range were created from the PIN data. The energy dependence of 
the pulse profiles in A 0535+26 are shown in Figure |2j The pulse profiles are complex in 
structure with narrow dips in the low energy ranges < 12 keV which morphed to become 
a simpler, more sinusoidal profile at higher energies. The following characteristics are 
observed with a careful examination of the profiles. 

1. A narrow dip at phase ~ 0.1, which decreases in strength with energy and disappears 
at energies > 14 keV. 

2. Indication of another sharp dip at phases ~ 0.2-0.3, which is evident only at the 
lowest energy range (< 2keV). 

3. The emission component between phases 0.5-0.7 becomes weaker and weaker with 
energy and finally disappears at ~ 17 keV. As a result the main dip of the profile (at 
phase ~ 0.6) is narrower at lower energies (< 12 keV) and broader at higher energies 
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'he energy 



f lNaik et al. 



depen dence is very similar to that found during the 2005 Suzaku observation 



20081 ). The profile is however, very different from th e simple sinusoidal 



at all energies found d uring the quiescence phase of the source (IMukherjee fc Pau 



profil e 



Negueruela et al. 



its giant outbursts ( IMihara 



2005 : 



2000) , or the double peaked profit 



1995 



Kretschmar et al 



e exte nding upto higher energies during 



19961) 



The energy dependence of the pulse profiles of XTE J1946+274 is shown in Figure |3j 
The pulse profiles show a clear double peaked structure which extends upto the high 
energies. The following characteristics can be observed in more detail. 

1. At the lowest energy ranges (0.3-4 keV), the peak (phase ~ 0.5) increases in strength 
with energy and the dip at phase ~ 0.8 increases in strength. 

2. Between 4-7 keV, the same dip mentioned above decreases in strength and the two 
peaks are almost equal in strength. 

3. Between 7-17 keV, this dip (phase ~ 0.8) disappears and a new, much weaker dip 
appears at ~ 0.9 which is probably the true interpulse region between the pulses. 

4. At the highest energies (25-70 keV), the second peak at phase ~ 0.1 becomes much 
weaker. 



The energy de pendence of the puls e profiles of XTE J1946+274 is very similar to that 



investigated by 



Wilson et al. 



(120031 ) during the 1998 outburst of the source. 
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3.3. Spectroscopy 



3.3.1. Pulse phase averaged spectroscopy 



We performed pulse phase averaged spectral analysis of A 0535+26 and XTE 
J1946+274 using spectra from the three front illuminated CCDs (XISs-0 and 3), the back 
illuminated CCD (XIS-1) and the PIN. We performed spectral fitting using XSPEC vl2.7.0. 
The XIS spectra were fitted from 0.8-10 keV and the PIN spectrum from 10-70 keV. 
The energy range of 1.75-2.23 keV was neglected due to an artificial structure in the XIS 
spectra around the Si edge and Au edge. After appropriate background subtraction, the 
spectra were fitted simultaneously with all parameters tied, except the relative instrument 
normalizations which were kept free. The XIS spectra were rebinned by a factor of 6 from 
0.8-6 keV and 7-10 keV, and by a factor of 2 between 6-7 keV. The PIN spectrum of A 
0535+26 was rebinned by a factor of 2 upto 22 keV, by 4 upto 45 keV, and 6 upto 70 keV. 
Due to comparatively inferior statistics in the PIN spectrum of XTE J1946+274, higher 
rebinning factors of 2, 6, and 10 were applied in the above mentioned energy ranges. 
In HMXB accretion powered pulsars, the continuum emission can be interpreted to 
arise by Comptonization of soft X-rays in the plasma above the neutron star surface. 



It is usually mode l ed phenomeno 



( IWhite et al. 



1983 



Mihara 



1995 



ogically with a powerlaw and cutoff at high energies 



Coburn 



200l|). The most widely used 



are the high energy cutoff (highecut) or the Fermi Dirac cutoff (fdcut) (JTanakalll986l ) with 



empirical mo dels 



the powerlaw component, or cutoff powerlaw (cutoffpl) mode 
negative-positive exponential powerlaw co mponent (NPEX) ( iMihara 



physical comptonization model 'CompTT' ( jTitarchuk 



. Other models include the 



19951 ). and a more 



1994J). We tried to fit the energy 



spectra with all the continuum models mentioned above, available as a standard or local 
package in XSPEC and carried out further analysis with only the models which gave best 
fits for the respective sources. 
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A 0535+26 : 



For A0535+26 the best fits were obtained with the NPEX, powerlaw and the 'CompTT' 
model (assuming spherical geometry for the comptonizing region). The powerlaw model 
however did not require a 'highecut' to fit the energy spectra. Including the GSO spectra 
in the fitting, the relative normalization of the GSO with respect to XIS showed that the 
flux in the GSO band (50-200 keV) was overestimated ~ 4 times without the inclusion 
of a 'highecut' in the spectrum. As inclusion of the GSO spectrum is not possible for 
phase resolved studies due to its limited statistics, and a spectrum of an accretion powered 
pulsar without a cutoff at higher energies is not viable, we have carried out further analysis 
with the 'NPEX' and 'CompTT' models. We applied a partial covering absorption model 
'pcfabs' in both the cases along with the Galactic line of sight absorption, to take into 
account the intrinsic absorption evident at certain pulse phases. This is evident in the pulse 
profiles and is a feature local to the neutron star. The narrow Fe ka feature found at 6.4 
keV was modeled by a gaussian line. In addition, a deep and wide feature found at ~ 45 
keV wa s modeled with a Lor entzian profile, which is the CRSF found previously in this 



source rtCaballero et al. 



20071 ). For A 0535+ 26, the CRSF has b een reported before at the 



same energy, even in a Suzaku observation (JTerada et al 



20061 ). So we do not comment 



on its detection significance here. We also tried a gaussian profile to model the CRSF 
feature. Since the centroid e nergy of the Lorentzia n description is not coincident with the 



minimum of the line profile (INakajima et al. 



20101 ). apart from a slight offset between the 



centroid energies of the Lorentzian and gaussian profile, the other parameters like the depth 
and width are consistent between the two models. The fits are also similar. We however 
considered a Lorentzian profile for the CRSFs for the rest of the paper after verifying the 
consistency between the Lorentzian and Gaussian profiles. The CRSF parameters were also 
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consistant within error bars for both the continuum models, the centroid energy being only 
slightly higher for the 'powerlaw' model. The reduced \ 2 obtained for the models were 1.25 
and 1.26 for 839 and 840 d.o.f respectively with no systematic residual pattern. 

XTE J1946+274: 

For XTE J1946+274, best fits with similar values of reduced x 2 were obtained with the 
'highecut', 'NPEX' and 'CompTT' model. Similar to A 0535+26, the local absorption of 
the neutron star was taken into account by the model 'pcfabs', and a gaussian line was also 
used to account for the narrow Fe ka feature found at 6.4 keV. A dee p and wide residua. 



Heindl et al. 



J2001I ). 



was found at ~ 38 keV, at the same energy as the CRSF discovered by 
As discussed previously, the CRSF was modeled with a Lorentzian profile. The 'highecut' 
and 'NPEX' models gave consistant values of the CRSF parameters, but the 'CompTT' 
model required a much shallower and narrow profile. Moreover, we were unable to constrain 
all the parameters of the 'CompTT' well for this source, probably due to the poorer quality 
of the PIN data. We have thus carried out the further analysis of this source with the two 
former models. For the best fitting models, the reduced \ 2 was 1-09 an d 111 respectively 
for 826 d.o.f. Without the inclusion of the CRSF, the difference in \ 2 was 150 and 119 
respectively for the same models. The best-fitting values for the spectral models for both 
the sources are given in Table 1. Figure [5] shows the best-fit spectra for both the sources 
along with the residuals before and after including the CRSF, thus showing the presence of 
the feature clearly. 



Miiller et al. 



(120121 ) however have reported the analysis of the RXTE, INTEGRAL and 
Swift observations during the same outburst of this source. Instead of a line at 36 keV, 
they found a weak evidence of a CRSF at ~ 25 keV. It may be worthwhile mentioning in 
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this context that the Suzaku PIN data has better sensitivity than INTEGRAL ISGRI at 
this energy range, and hence may be better suited for CRSF detection. However we have 
carefully checked the statistical significance and possible systematic errors associated the 
CRSF. 

Statistical significance: To estimate the detection significance of the CRSF we tried 
to fit the PIN spectrum alone with the 'highecut' model with its powerlaw index frozen to 
the value obtained from the best fitting broadband spectrum. The addition of the CRSF 
improved the \ 2 from 51.56 to 28.13 for 20 d.o.f corresponding to an F value of 16.7, and a 
F-test false alarm probability of 6 x 10~ 4 . 

Possible systematic errors: At first, we used the the earth occultation data to check 



the reproducibility of the NXB ( iFukazawa et al. 



20091 ). We extracted the spectra using the 



earth occultation data in three energy bands centering the CRSF and compared ratio of the 
count rates with the NXB. The ratio obtained were 1.3, 1.2 and 1.2 at 10 — 28, 28 — 48 and 
48 — 70 keV respectively indicating the lack of any energy dependent feature that can be 
introduced by the simulated X-ray background. We also included a systematic uncertainty 
of 3% on the PIN spectrum to check the detection of the CRSF. The line was still detected, 
but the uncertainty in the depth of the feature increased by 23%. The detection of pulse 
phase dependence of this feature as discussed in section 3.3.2 is also in favor of its presence 
since the background data is not expected to vary over the pulse phase. Finally, to verify 
the existence of the CRSF in a model independent manner, we divided the PIN spectrum 
of a pulse phase with the deepest CRSF, by the same of a pulse phase with the shallowest 
CRSF detected (see section 3.3.2, pulse phase resolved spectroscopy for the corresponding 
spectra). Figure 0] shows the ratio plot of the two spectra. Although the quality of the data 
is not good after 40 keV, the dip at ~ 30-35 keV is clearly seen indicating the presence of 
the CRSF. 
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3.3.2. Pulse phase resolved spectroscopy 

For the phase-resolved analysis we extracted the source spectra for both the XIS's and 
the PIN data after applying phase filtering in the FTOOLS task XSELECT. The same back- 
ground spectra and response matrices as used for the phase- averaged spectra were however 
used in both the cases. The spectra were also fitted in the same energy range and rebinned 
by the same factor as in phase- averaged case. The Galactic absorption (A/hi) column density 
and the Fe line width were frozen to the phase-averaged values for the two respective models. 



Phase resolved spectroscopy of the cyclotron parameters : For investigating 
the pulse phase-resolved spectroscopy of the two CRSFs, phase resolved spectra were 
generated with their phases centered around 25 independent bins but at thrice their widths. 
This resulted in 25 overlapping bins out of which only 8 were independent. We however 
froze the width of the CRSF to the phase-averaged value of the respective models, and 
varied the rest of the continuum as well as the line parameters with pulse phase. This was 
due to our inability to constrain all the parameters because of limited statistics. Figure |6] 
shows the variation of the cyclotron parameters of the sources using the best fit models as 
a function of pulse phase. For both the sources, the different continuum models used result 
in a very similar pattern of variation of the parameters. This gives us a reasonable amount 
of confidence on the obtained results. The following features are evident from the Figure |6j 
The results are compared with respect to the high energy PIN profile (10-70 keV). 
A 0535+26 : 

1. the energy (El cyc i) varies by 14% (~ 43-50 keV). The pattern of variation of both 
the energy (£l cyc i) and depth (Dl cyc \) has a gradually increasing trend with the pulse 
profile and drops off abruptly in the off-pulse region (phase ~ 0.6), picking up again 
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where the pulse profile picks up. 

2. The depth (El cyc \) cannot be constrained at all phases by both the models, and at 
the off pulse phase at ~ 0.6, only the 'CompTT is able to constrain the depth. It has 
a very sharp pattern of variation, varying between ~ 0.8-4, and it is shallowest near 
the pulse peak and deepest near the pulse minima. 

XTE J1946+274 

1. The energy (El cyc i) varies about 36%. It's value is generally higher in the first pulse 
with the values peaking near the first peak (phase~ 0.7-0.8), and a decreasing trend 
near the second pulse. 

2. The depth (ET cyc i) varies between 1-3. It is deepest at the interpulse regions at phase 
~ 1.0 and shallow between phase 0.5-0.8 near the first peak. Due to limited statistics, 
specially of the PIN spectra, the CRSF parameters however cannot be constrained at 
the main dip, and at the ascending edges of the first peak (phase ~ 0.5-0.7). 

Phase resolved spectroscopy of the continuum parameters: A dependence 
of the continuum energy spectrum on the pulse phase is implied from the strong energy 
dependence of the pulse profiles, as seen in Figure [2] and Figure [31 A partial covering 
absorption model in which the absorber is phase locked with the neutron star is required 
to explain the narrow energy dependent dips in the pulse profiles. This was also our main 
motivation in applying the partial covering absorption 'pcfabs' to model the continuum 
energy spectra. We generated the phase resolved spectra with 25 independent phase bins 
to investigate the pulse phase-resolved spectroscopy of the continuum parameters for A 
0535+26. Due to the short spin period of XTE J1946+274, 25 independent phase bin 
extraction was not possible, specially for the XIS data. We proceeded with the extracting 
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of 25 overlapping but 8 independent phase bins for extraction of both XIS and PIN data 
as was done for the phase resolved spectroscopy of the CRSF parameters. The cyclotron 
parameters of the corresponding phase bins were frozen to the best-fit values obtained from 
the results of investigation of the cyclotron line parameters using 25 overlapping phase 
bins. Figures |8] & [9] shows phase resolved continuum parameters using the best-fit spectral 
models as a function of the pulse phase for A 0535+26 and XTE J1946+274 respectively. 
The results obtained as seen from the Figure from both the models are as follows: 

1. For both the sources, there is an abrupt increase in the value of the local absorption 
component (iV^), with a corresponding change in the value of the covering fraction 
(Cvf rac t) at the dips of the low energy XIS profile. This picture is in agreement with 
a narrow stream of matter present at those phases responsible for absorption of the 
low energy photons. The properties of the plasma in the accretion stream, which may 
be a narrow structure having different values of opacities and optical depths can be 
traced from the changes in the value of A^H2 and the covering fraction. As can also 
be seen clearly, the main strength of our results lie in the fact that we have obtained 
similar patterns of variation of Nn 2 and Cff ract using different continuum spectral 
models for the sources. 

2. There are also corresponding changes in the other continuum parameters like the 
powerlaw photon index (r) of the 'powerlaw', seed temperature (CompTTi ), optical 
depth (r) and KT of the 'CompTT' model for A 0535+26, and the powerlaw photon 
index (r), the E-folding and E-cut energy of the 'highecut' and the NPEX al and 
'KT' of the 'NPEX' model for XTE J1946+274. The main aim of this paper is 
however the pulse phase resolved variation of the CRSF parameters, and detailed 
discussion of these results are beyond the scope of this paper. 
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4. 4U 1907+09 

4.1. Timing analysis: Light curves & Hardness ratio & pulse period 

determination 



4U 1907+09 is a variable X-ray source showing flaring and dipping activity in the 



timescales of minutes to hours (lin 't Zand et al 



19971 ). We performed timing analysis after 



applying barycentric corrections to the event data files using the FTOOLS task 'aebarycen'. 
Light curves were extracted with a time resolution of 8 s (full window mode of the XIS 
data) for the XISs (0.2-12 keV), and 1 s for the HXD/PIN (10-70 keV) respectively. We 
applied pulse folding and y 2 maximization technique to search for pulsations in the XIS 
data. The source having an eccentric orbit with a short orbital period, proper correction of 

se period. However, t he 



in 't Zand et al. 



19981 ) . 



the pulse arrival times are required to accurately determine the pu. 

orbital ephemeris of this source is not known with high accuracy (ji 

Thus to account for the orbital motion of the binary, we included a ■£ term in the fitting, 

starting with an initial guess consistant with the parameters of the binary, and iterating for 

different values of -£ to get the maximum % 2 . The best fit period corresponding to this was 

441.113 ±0.035 s MJD 542 09.43189 with jg = 3.1 x 10~ 6 . This value obtained is marginally 



higher than that found by 



Rivers et al 



<|2O10h (441.03 ± 0.03). However they have not 



mentioned, taking into account the orbital correction of the pulse arrival times in their work 
which might be a reason for this discrepancy. Figures [1] shows the XIS and PIN light curves 
along with the hardness ratio. As can be seen from the figure, the light curves show two 
flaring features in betw een and a dip in the last ~ 10 ks of the observation. These features 



were also mentioned in iRivers et al. 



(120101 ). while performing time resolved spectroscopy 
of the same Suzaku observation, and were probed further by them to investigate the 
spectral variability with time. The flares may, however also affect our results of pulse phase 
resolved spectroscopy. We have thus compared the pulse profiles and the energy spectra in 
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these stretches individually with that from the rest of the observation. Though the pulse 
profiles look very similar in all the stretches, the energy spectra is harder with an increased 
absorption in the last stretch of the observation containing the dip. The main aim of this 
work being pulse phase resolved spectroscopy to probe the CRSF parameters, we excluded 
the stretch of the observation coincident with the dip in the light curve for further analysis. 
The arrows in Figures [T] indicate the length of the observation chosen for this work. Pulse 
profile for this duration of observation was also created in the XIS and PIN energy bands as 
before for A 0535+26 and XTE J1946+274. Due to the absence of low energy dips in this 
source however, the energy dependence of the pulse profiles was not investigated further. 

4.2. Pulse phase averaged spectroscopy 



Phase averaged spectroscopy was carried out in the same procedure as in A 0535+26 
and XTE J1946+274. Best fits were obtained with the 'highecut', 'NPEX' and 'compT 



mode 



with comparable values of reduced y 2 and similar residual patterns. 



Rivers et al. 



(120101 ) also obtained similar results with the 'highecut', 'fdcut' and 'NPEX' model. A 



comparison between the NPEX model parameters obtained in our analysis and those 



reported in 



Rivers et al. 



( 120101 ) reveal a softer less absorbed spectra obtained by us. This 
is expected, since we have excluded the the last stretch of data from our analysis which 
had a more harder and absorbed spectra. Two gaussian lines were also used to model 
the narrow Fe ka and Fe k/3 feature found at 6.4 and 7.1 keV respectively. In addition, a 
relatively shallow and narrow feature found at ~ 18 keV was modeled with a Lorentzian 



profile which is the CR SF previously detecte d in this sour c e (|Makishima fe Mihara 



Makishima et al. 



1999). As also discussed in 



Rivers et al. 



1992 



(120101 ). the first harmonic of 



the CRSF at ~ 36 keV could not be detected in the PIN spectra probably due to the 
statistical limitation of the data in this energy range. The CRSF parameters obtained with 
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the 'NPEX' and 'C ompTT' models were consistant within error bars with that found by 



Rivers et al. 



(120101 ) who performed phase resolved spectroscopy in 6 independent bins using 



the gaussian absorption model (keeping in mind that the centroid energy of the Lorentzi an 



20101)). 



description is not coincident with the minimum of the line profile ( iNakajima et al. 
The 'highecut' model however required a deeper CRSF to fit the spectra. The reduced \ 2 
obtained for the models were 1.62, 1.51 and 1.69 for 832, 837 and 838 d.o.f for highecut, 
NPEX and CompTT respectively with no systematic residual patterns. Due to the 
compatibility of the CRSF parameters obtained with the 'NPEX' and CompTT models, we 
have carried out further phase resolved analysis using these two models. Figure [5] shows 
the best-fit spectra for 4U 1907+09 along with the residuals before and after including the 
CRSF, thus showing the presence of the feature clearly. The CRSF in this source is ver y 



stron g and has also been reported in the same Suzaku observation before (iRivers et al. 



2010). We therefore do not comment on its detection significance. 



4.3. Pulse phase resolved spectroscopy 

For investigating the pulse phase-resolved spectroscopy of the CRSF, we generated 
phase resolved spectra with 25 overlapping but 8 independent phase bins and used the 
same analysis procedure as discussed previously for A 0535+26 and XTE J1946+274. We 
were however able to constrain the phase dependent variation of all the CRSF parameters 
for this source, probably due to the longest observation duration available for it and a 
low cyclotron energy compared to the other sources. Figure [6] shows the variation of the 
cyclotron parameters of the source using the best fit models as a function of pulse phase. 
As in the case of the previous sources, the similar pattern of variation obtained for the 
different continuum models used give us considerable amount of confidence on the obtained 
results. The following characteristics can be observed in more detail from Figure [61 As 
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before, the variations are compared with respect to the high energy PIN profile. 

1. The energy El cyc i varies by ~ 19%. Its value is maximum near the peak of the first 
pulse (20 keV at phase ~ 0.3), and again at the ascending edge of the second pulse 
(phase ~ 0.6), the minimum being at the second pulse peak (15 keV at ~ phase 
0.7-0.8). 

2. The depth -Elcyci has a clear double peaked pattern with the peaks corresponding 
to the ascending edge of the first pulse and the peak of the second pulse (phase 

~ 0.1 — 0.2 and 0.7 respectively). It is minimum near the pulse minima (phase ~ 0.9). 
-El cyc i varies between 0.2-1.4 and is generally greater for the first pulse. 

3. The width (Wl cyc \) has a similar pattern of variation as £T cyc i, and peaks at similar 
phases with values varying within 7 keV. 

5. Discussions & Conclusions 

In the present work we have presented the results of detailed pulse phase resolved 
spectroscopy of the CRSF parameters of A 0535+26, XTE J1946+274 and 4U 1907+09 
using long Suzaku observations. Pulse phase dependence of the CRSF parameters are 
obtained for the first time in A0535+26 and XTE J1946+274 and a more detailed and 
careful analysis has been done in 4U 1907+09 which is consistant with the earlier result 



obtained using the same observation (IRivers et al. 



2010j). The analysis is done taking into 



account various factors which might smear the pulse phase dependence results as mentioned 
in earlier sections. The strength of our results lie in the fact that we have obtained similar 
pattern of variation of the CRSF parameters for all the three sources with more than one 
continuum model. 
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5.1. pulse phase dependence of the cyclotron parameters 



Results of pulse phase resolved spectroscopy of the cycl otron parameters have been 



presented previously in som e sources, for e xample in Her X-l (ISoong et al 



2008 



1999 



Kloch 



2002 



•cov et al. 



2008), 4U 0115+63 ( 



La Barbera et al. 



Cen X-3 (ISuchv et al. 



( jSuchy et al. 



2011; 



2003 



Heindl et al 



Maitra fc Pau 



)0) 



2000), Vela X-l 



2013|), 4 U 1538-52 ( 



2008 ), and more recently in GX 3 01-2 (ISuchv et al. 



Maitra et al. 



1990 



Enoto et al. 



Krcvkcubo hm et al. 



Robba et al. 



20121 ) and 4U 1626-67 fllwakiri et al. 



2012)). 



20011), 



20121 ). 1A 1118-61 



By modeling the pulse phase dependence with different continuum models, we have been 
able to establish the robustness of the results. In the process of trying to fit the energy 
spectrum with different continuum models we have also noticed certain trends in continuum 
model fitting. The 'highecut' model being a very simple model with less number of 
parameters, is a good choice to model the continuum in case of moderate or poor statistics. 
This is evident in the case of XTE J1946+274. The 'CompTT' on the hand which is a more 
physical description of the spectra and has a reasonable number of free parameters is better 
for continuum fitting specially for phase resolved spectroscopy if the statistical quality of 
the data is reasonably good. This is probably the reason why it failed to constrain the 



continuum well in the case of XTE J1946+274. The 'NPEX' model approximates 



the 



photon number spectrum for an unsaturated Comptonization (jSunyaev fc Titarchuk 



Mcszaros 



19801 : 



1992J), and has a clear physical meaning inspite of being a phenomenological 



model. It is useful for all the three sources with significantly different statistical quality. 



By assuming certain physics and geometry of the line for ming region, the CRSF 



(120081 . 



2Qll|); 



Arava fc Hardin g (1999); 



feature has been modeled analy ti cally and with simulat ions by 

^v ^, 

Arava-Gochez fc Harding (2000); Schonh err et al.l (120071 ) and more recently by 



Nishimura 



Mukherjee fc Bhattacharyal (120121 ). Although these models predict variations 



in the depth, width and the centroid energy of the CRSF features with the changing 
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viewing angle at different pulse phases, a variation in the CRSF parameters as large as 
30% as found in our results needs to take into account either a possible deviation or 
distortion from the simple 



dipole geometry of the magnetic field (Schon herr et al 



Mukherjee fc Bhattacharya 



20121 ) , a gradie nt in t 



different geometry of the accretion column (JKraus 



re field itself (INishimura 



2007; 



200$), or a 



200 ll ). A detailed modeling taking into 



account these factors would provide us a detailed information about the geometry and 
emission patterns of the sources. However simpler interpretations can be done, since the 
correlation of the deepest and shallowest CRSFs with the pulse profile of the source can 
provide some idea about the beaming pattern of the source at that luminosity. Following 
this, the trend of shallowest lines near the pulse peak and deepest near the off-pulse as 
found in A 0535+26 and XTE J1946+274, favors a pencil beam geometry. On the other 
hand, deepest and widest lines found near the peak and shallowest and narrowest near 
the off-pulse as found in 4U 1907+09 favors a fan beam geometry for the emission. These 
results may be further complicated by assuming the contribution from both the magnetic 
poles of the neutron star in contrary to one of them, either due to gravitational light 
bending or particular geometry of the system allowing the view of both the poles. Modeling 
of the variations of the CRSF parameters with pulse phase is ongoing. Detailed discussions 
on the same will be made in a future work. 



This research has made use of data obtained through the High Energy Astrophysics 
Science Archive Research Center On line Service, provided by NASA/Goddard Space Flight 
Center. Chandreyee Maitra would like to thank Carlo Ferrigno for providing the 'fdcut' 
and the 'newhighecut' local models. 
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Fig. 1.— Light curves of A 0535+26, XTE J 1946+274 and 4U 1907+09 obtained with 
Suzaku.The first panels in each figure shows the light curve for one of the XIS in the energy 
band of 0.3-12 keV. The second panel shows the same obtained in the PIN energy band(10- 
70 keV). The time binning is equal to the respective pulse periods for A0535+26 and 4U 
1907+09 and 10 pulsar period in the case of XTE J1946+274. The bottom panel shows 
the hardness ratio. The arrows in the hardness ratio of 4U 1907+09 indicate the stretch for 
which data was chosen to perform phase resolved analysis. 
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Fig. 2. — Energy dependent pulse profiles of A 0535+26 using XIS & PIN data. The energy 
range for the pulse profiles are specified inside the panels. 
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Fig. 3. — Energy dependent pulse profiles of XTE J1946+274 using XIS & PIN data. The 
energy range for the pulse profiles are specified inside the panels. 
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Fig. 4. — Ratio of counts of the energy spectrum with the shallowest cyclotron line to counts 
of the spectrum with the deepest cyclotron line in XTE J1946+274. Though the ratio of the 
counts after 40 keV have large error bars due to statistical limitations, the dip in counts at 
~ 30-40 keV is clealy visible indicating the presence of the CRSF 
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Fig. 5. — The pulse-phase averaged spectrum of A 0535+26, XTE J1946+274 and 4U 
1907+09 showing all the individual model components (Top left, right and bottom respec- 
tively). The upper panel shows the best-fit spectra as obtained with the 'NPEX' model. 
The middle panel shows the residuals without inclusion of Lorentzian profile for the CRSF 
in the spectra, and the bottom panel shows the residuals after the inclusion of the CRSF. 
This clearly shows the presence of the CRSFs in the energy spectra of all the sources 



-34- 



1.5 



- 1 



0.5 



Q. 



2 



o.x 



50 






40 
6 




I I I I 



f-H — I — I — | — I — I — I — I — |- 



,+ 1 



■l-H — I — I — \ 



I ' ' ' ' I ' 



IE 



i— | — l-H — l-H — |— f 



W#*#W^M 



m 



H 



1.4 
1.2 



Oh 



on 

" 0.8 

0.6 

1.2 



Oh 1 

Oh 



1 - 



Q.8 



0.6 
50 



45 



~H> 

WJ^40 



35 

4 
3 

2 
1 





■ I TH i J 1 1 i i i i | Tm i i' ' | i i i i i 





0.5 1 1.5 



Pulse phase 



0.5 1 1.5 2 

Pulse phase 

Fig. 6. — Variation of the cyclotron line parameters in A 0535+26 (left panel) and XTE 
J1946+274 (right panel) as is obtained with the two models. In the left panel, the black 
points denotes the parameters as obtained with the 'NPEX' model The red points denote 
the parameters as obtained with the 'CompTT' model. In the right panel the black points 
are obtained with the 'highecut' model and the red points with the 'NPEX' model. Only 8 
of the 25 bins are independent. The XIS (0.3-10 keV) and PIN (10-70 keV) pulse profiles 
are shown in the top two panels respectively which denote the normalized intensity. 
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Fig. 7. — Variation of the cyclotron line parameters in 4U 1907+09 as is obtained with the 
two models. The black points denotes the parameters as obtained with the 'NPEX' model 
The red points denote the parameters as obtained with the 'CompTT' model. Only 8 of the 
25 bins are independent. The XIS (0.3-10 keV) and PIN (10-70 keV)pulse profiles are shown 
in the top two panels respectively which denote the normalized intensity. 
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Fig. 8. — Variation of the spectral parameters with phase along with the pulse profile (0.3-12 
keV for XIS and 10-70 keV for PIN) in A 0535+26. The left panel shows the variation using 
the model 'NPEX' and the right panel shows the same using the model 'CompTT'. The pulse 
profiles denote the normalized intensity. The XIS (0.3-10 keV) and PIN (10-70 keV) pulse 
profiles are shown in the top two panels respectively which denote the normalized intensity. 
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Fig. 9. — Variation of the spectral parameters with phase along with the pulse profile (0.3-12 
keV for XIS and 10-70 keV for PIN) in XTE J1946+274. The left panel shows the variation 
using the model 'highecut' and the right panel shows the same using the model 'NPEX'.The 
pulse profiles denote the normalized intensity. The XIS (0.3-10 keV) and PIN (10-70 keV) 
pulse profiles are shown in the top two panels respectively which denote the normalized 
intensity. 
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Table 1: Best fitting phase averaged spectral parameters of A 0535+26, XTE J1946+274 
and 4U 1907+09. Errors quoted are for 99 per cent confidence range. 



Model 


A 0535+26 


XTE J1946+274 


4U 1907+09 


parameters 


NPEX 


CompTT 


highecut 


NPEX 


NPEX CompTT 


iVg-L (10 22 atoms cm,-' 2 ) 


0.59 ±0.02 


S1+ ' 04 


1.30 ±0.028 


1 9 7 +0.03 
1 - z '-0.04 


1.97 ±0.01 1.62 ±0.03 


N^ 2 (10 22 atoms cm,- 2 ) 


q 7*1+0.25 
°- ,y: -0.36 


741 +1.45 


5-86_ ; 82 


9 ns+ 2 - 72 


- - 


C^fract 


0.49 ±0.01 


0.37 ±0.02 


0.27 ±0.04 


17+0- 04 


- - 


Powlndex 


- 


- 


1 09+0- 05 


- - 


E-folding energy (keV) 


- 


- 


95 57+ 2 ' 75 - 
zo -°'-2.26 


- - 


E-cut energy (keV) 


- 


- 


7 n r,+0.69 
'■ uz -0.29 


- - 


powerlaw^ OTIn 


- 


- 


0.021 ± 0.001 


- - 


CompTT T0 (keV) 


0.52 ±0.04 


- 


- 


0.47 ±0.01 


CompTT KT (keV) 


14.67+.?;°? 


- 


- 


4.78 ±0.04 


CompTT t 


7 ,,+0.44 
'•J-J_0.48 


- 


- 


14.87 ±0.08 


CompTT nolrn c 


0.008 ± 0.0009 


- 


- 


0.0317 ± 0.0005 


NPEX al 


, nl +o.05 _ 


- 70+0- 08 

u -' u -o.io 


0.84 ±0.02 


NPEX a2 


-2.0 (frozen) 


-2.0 


(frozen) 


-2.0 (frozen) 


NPEX KT (keV) 


11.49ti;« 


13.2912:74 


4.41 ±0.04 


NPEX nolln l c 


0.041±°;°°? 


- 014+ 0002 


0.072 ±0.001 


NPEX nolm 2 C 


4.94e - 6 ± 1.25e - 6 


2.34 d ± 1.5 


2.47e - 4 ± 7.5e - 6 


^Icycl 


42 60+ ' 91 


4S 24+ ' 85 


38.30++^ 


38.65+1;% 


17.96+.Q; 20 , 18.07 ±0.18 


Slcycl 6 


i q 7 +0.25 
1 -°'-0.23 


I 4 q+0.35 
1 - 41o -0.26 


1 79 +0.41 
1 -' z -0.28 


1-50+clJ! 


0.68 ±0.03 0.61 ±0.03 


Wl cycl 


7 1 7+2.28 
'• 1 '-1.82 


4 Q3+ 1 - 37 


q fi i+3.69 


8 .6i±i:« 


, ,,,+0.39 9 7 o+0.39 
■J-^-O^ z -'°-0.37 


Iron line energy (keV) 


- 


- 


- 


- 


71n +0.035 71n +0.033 
'• iu -0.024 '- iu -0.037 


Iron line eqwidth (eV) 


- 


- 


- 


- 


10.74 ±2.90 7.9 ±2.6 


Iron line energy (keV) 


6.41± 0.01 


6.41± 0.01 


6.41 ±0.022 


6.41 ±0.02 


fi 42 +0.008 fi 42 +0.007 
"•* z -0.010 D - 4z -0.009 


Iron line eqwidth (eV) 


23.37 ±4.61 


24.28 ±5.26 


29.05 ±4.93 


28.39 ±5.15 


51.67 ±3.62 43.81 ±1.66 


Flux (XIS) f (0.3-10 keV) 


3.05 ±0.04 


3.05 ±0.05 


1.80 ±0.02 


1.80 ±0.02 


4.75 ±0.06 4.74 ±0.06 


Flux (PIN) 9 (10-70 keV) 


8.14 ±0.02 


8.14 ±0.03 


3.47 ±0.01 


3.59 ±0.0 


6.87 ±0.02 6.91 ±0.02 


reduced x 2 /d.o.f 


1.25/839 


1.26/837 


1.09/826 


1.11/826 


1.51/837 1.69/838 



a Denotes the Galactic line of sight absorption b Denotes the local absorption by the partial covering absorber 'pcfabs' 
c photons keV -1 cm~ 2 s _1 at 1 keV 



*10" 



; c denotes the optical depth of the feature. 



10 ° ergs cm 2 s x and are in 99 % confidence range. 



' 10 10 ergs cm 2 i 



and arc in 99 % confidence 



